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Abstract 

We incorporate chiral symmetry breaking in a soft- wall version of the AdS/QCD 
model by using a modified dilaton profile and a quartic term in the bulk scalar 
potential. This allows one to separate the dependence on spontaneous and 
explicit chiral symmetry breaking. The resulting mass spectra in the scalar, 
vector and axial-vector sectors compares favorably with the respective QCD 
resonances. 
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1 Introduction 

The AdS/CFT correspondence [1, 2, 3, 4] provides a remarkable mathematical tool 
with which to understand strongly-coupled gauge theories. Calculations performed 
on the gravity side can be reinterpreted as due to nonperturbative effects on the 
field theory side. The theory of the strong interactions, quantum chromodynamics 
(QCD), is a strongly-coupled gauge theory at low energies, and therefore a natural 
candidate with which to apply the gauge/gravity correspondence. This had led to 
a bottom-up approach, commonly known as AdS/QCD [5, 6], which relates QCD 
to a five-dimensional (5D) gravity theory. This simple model captures the essential 
features of the low-lying meson spectrum. 



To incorporate more realistic features of the excited states, the AdS/QCD model 
can be modified to include a dilaton with a quadratic profile [7]. While the linear 
radial spectrum is indeed produced in this soft- wall version of the AdS/QCD model, 
the form of chiral symmetry breaking is not QCD-like. In particular, the bulk scalar 
field, dual to the quark bilinear operator qq, whose vacuum expectation value is 
responsible for spontaneous chiral symmetry breaking, does not allow the spontaneous 
and explicit breaking to be independent. Moreover, chiral symmetry is restored for 
the highly excited states, a phenomenological feature that is not supported in the 
QCD mass spectrum [8]. 

We [9] modify the existing soft-wall version of the AdS/QCD model in order to 
incorporate these two phenomenological features of QCD. This is done by adding 
a quartic term to the bulk scalar potential and changing the dilaton profile. By 
assuming that the bulk scalar field contains the desired limiting behavior for non- 
restoration of chiral symmetry [4, 8], we derive a new dilaton background profile. 
The extra parameter introduced by the quartic term decouples the quark mass from 
the chiral condensate, thereby allowing for spontaneous and explicit chiral symmetry 
breaking to occur independently. 



2 The Model 



We consider a modified version of the soft-wall AdS/QCD model first introduced in 
Ref. [7] and investigated further in Refs. [10, 11, 12, 13, 14, 15, 16]. The background 
geometry is assumed to be 5D AdS space with the metric 



ds 2 = gMN<ix M dx N = a 2 (z) {r] liv dx^dx l/ + dz 2 ^ 



where a(z) = L/z is the warp factor, L is the AdS curvature radius and the Minkowski 
metric 7]^ = diag(— 1, +1, +1, +1). The conformal coordinate z has a range < z < 
oo. To obtain linear trajectories, Ref. [7] also introduced a background dilaton field 
4>(z) = \z 2 where A sets the mass scale for the meson spectrum. 

To describe chiral symmetry breaking in the meson sector the 5D action contains 
SU(2) L x SU(2)i? gauge fields and a bifundamental scalar field X. As suggested by 
Ref. [7], we add a quartic term in the potential V(X) of our 5D action, 
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where m 2 x = —3/L 2 , k is a constant, g\ = Yin 2 /JV C , and the covariant derivative is 
D M X = d M X - iA^X + iXA^. To describe the vector and axial-vector mesons we 
simply transform to the vector (V) and axial- vector (A) fields where V M = \{A^ + 
A%) and A M = |«-<). 

The scalar field X, which is dual to the operator qq, is assumed to obtain a 
z-dependent vacuum expectation value (VEV), 
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which breaks the chiral symmetry SU(2) L xSU(2) R toSU(2)y. Assuming (3) we obtain 
a nonlinear equation for v(z), 



d*(a 3 e _ *d*u(*)) -a h e-' l> {m 2 x v{z) 



-v\z)) = 0. 
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When k — 0, the solution of (4) which leads to a finite action in the limit z — > oo is 
given by [7, 14] 



where U (a, b, y) is the Tricomi confluent hypergeometric function. 

As expected from the AdS/CFT dictionary established in Refs. [4, 3], as z — > 
the VEV should take the form 



Here m q is the quark mass, S = (qq) is the chiral condensate, and ( is the normaliza- 
tion parameter introduced in Ref. [13]. For fixed values of m q and X, the introduction 
of ( still satisfies the Gell-Mann-Oakes-Renner relation, m^f% = 2m 9 S. In actuality 
the normalization ( is not a free parameter but is determined by QCD, as shown in 
Ref. [13], to be C = V / 3/fi , 5- Expanding the solution (5) in the small z limit leads to 
a oc m q and (3 oc £ oc m q . Thus in the limit m q — > the model eliminates explicit 
and spontaneous chiral symmetry breaking in contradiction with QCD. It will be seen 
that the introduction of a quartic term in the potential V(X) avoids the dependence 
of the chiral condensate on the quark mass encountered in Refs. [7, 12, 14]. 

The solution (5) has an asymptotic limit v(z) — > constant for large values of z. 
This asymptotic behavior suggests chiral symmetry restoration in the mass spectrum, 
a phenomenon not supported in QCD (although speculation continues on whether 
such a restoration indeed exists [17, 18]). As noted in Ref. [8] the highly excited 
mesons exhibit seemingly parallel slopes signifying that chiral symmetry is not re- 
stored. In order to incorporate this behavior v(z) must behave linearly as z becomes 
large, v(z — > oo) ~ z causing the mass-squared difference between vector and axial- 
vector resonances to approach a constant as z — > oo. By including a quartic term and 
requiring v(z) to have this linear asymptotic behavior we aim to incorporate these 
Q CD-like characteristics into the soft- wall model. 

Instead of solving for v(z) directly, we assume the VEV asymptotically behaves 
as expected, namely Eq. (6) and 
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and then solve for the dilaton <f>(z) using (4), which becomes 
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where the prime denotes the derivative with respect to z. Given the required behavior 
(6) and (7) we can uniquely determine the dilaton profile up to a constant. With this 



procedure the two sources of chiral symmetry breaking decouple while simultaneously 
allowing for linear trajectories in the meson spectrum. 

A particularly simple parametrized form for v(z) that satisfies (6) and (7) is 

v (z) = y(A + BtanhCz 2 ), (9) 

The parameters 7, A, B and C can be expressed in terms of the input parameters 
m q , E, A, k as 
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The parameter A is determined by the average slope of the radial trajectories of the 
scalar, vector, and axial- vector mesons for radial quantum numbers n > 3. Its value 
was determined to be A = 0.1831 GeV 2 , as explained in the next section. The quark 
mass, quark condensate, pion decay constant, and pion mass are all related by the 
Gell-Mann-Oakes-Renner relation, f%m% = 2m g E. This relation holds in this model 
as a natural consequence of chiral symmetry [5]. We use the measured values of 
f n = 92.4 MeV and m n = 139.6 MeV, and adjust the quark mass to reproduce the 
input value of f n from a solution to the axial-vector field equation in Section 4 for 
a given value of k. The parameter k essentially controls the mass splitting between 
the vector and axial- vector mesons. It is determined to be k — 15 by a best fit to 
the radial spectra of the axial- vector mesons. This results in m q = 9.73 MeV and 
therefore E = (204.5 MeV) 3 . The inferred value of the quark mass is consistent with 
an average of the up and down quark masses as summarized in the Review of Particle 
Physics [19] as appropriate at the hadronic energy scale. 

Other parameterizations of v(z) which lead to qualitatively similar behavior as 
that required in (6) and (7) can also be used. Other forms for v{z)/z include (01 + 
a2-2 2 )/(l + as? 2 ) [Pade] and c\ — C2 exp(— c%z 2 ) [Gaussian]. These forms were all 
studied but the best results were found to be obtained using the form (9). These 
parameterizations are shown in Fig. 1 using the above parameters. The resulting 
dilaton profile 4>(z) is shown in Fig. 2. 

3 Meson Mass Spectra 

The soft wall model can be used to fit the observed scalar, vector, and axial-vector 
meson mass spectra [19, 20, 21, 22]. All but one of the included states are listed 
in the Review of Particle Physics (RPP) [19]. For an explanation of the masses 
used see Ref. [9]. A straight line is fitted to the m 2 versus n plot with n > 3 




Figure 1: A plot of v(z)/z for various parameterizations fitted to the mass spectra. 
The best fit to the mass spectra is obtained with the tank form (9). 
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Figure 2: The dilaton profile (fi(z) resulting from the various parameterizations of 
v(z). The best fit to the meson spectra occurs with the tanh parameterization (9). 
For z < 1 the behavior deviates from the quadratic asymptotic form. 
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Figure 3: A straight-line fit to the measured scalar, vector and axial-vector mass 
spectra for n > 3 used to determine the dilaton mass parameter A. 



for all three mesons, assuming the same slope 4A but different intercepts, namely 
ml = 4An + m\. The results are: A = 0.1831 ± 0.0059 GeV 2 , m 2 Vfl = 0.0806 ± 0.0104 
GeV 2 , m\ = 1.5023 ± 0.0366 GeV 2 , and m 2 = -0.6634 ± 0.0038 GeV 2 . See Fig. 
3. We use this value of A in our model calculations. 

3.1 Scalars 

Introducing a quartic term in the Lagrangian causes the scalar excitations to couple 
with their own VEV, giving a modified equation of motion unlike those in Refs. [6, 14]. 
The equation of motion can be brought into a Schrodinger-like form 



where u s = <f>(z) + 3 log 2. To solve this equation we implement a shooting method in 
which (14) is solved for various values of m Sn . The eigenvalues are then those mass 
values that produce a solution for S n (z) that is bounded as z — > oo. Applying the 
shooting method to (14) with the boundary conditions lim Zo _> s n( z o) = 0, d z s n (z — > 
oo) = produces a scalar mass spectrum which has the correct slope but too large 
magnitude. This could well be a failure of this specific model. However, it has been 
argued [23] that the a meson is not truly a particle but only an enhancement in the 
7r-7r cross section, and that the phase shift does not reach 90°. Removing the / (550) 
would shift all the higher masses to the left by one unit of n, resulting in a better 
fit to the model. See Fig. 4. An obvious extension of this work would be to include 
strange quarks and glueballs and to determine the mixing among the resulting scalar 
states. 
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Figure 4: Comparison of the predicted scalar mass eigenvalues using the tanh form 
(9) ofv(z) (solid) with the QCD f scalar mass spectrum with the a meson removed 



From the action (2) the equation of motion of the vector field can be written in the 
Schrodinger form 



where ui = (f>(z) + log z. Using the dilaton form = \z 2 , the eigenvalues of (15) at 
large n can be solved analytically with the boundary conditions \im Z0 ^ Q v n (z ) = 0, 
d z v n (z — > oo) = and agree with those found in Ref. [7]. However, since the dilaton 
specified in (8) is modified for z < 1 there is a change in the slope of the mass 
spectrum around n = 2 which matches the behavior of the experimental data. The 
numerical vector mass spectrum is compared to the experimental data in Fig. 5. 
While the prediction for the p(775) mass is low, the rest of the vector meson masses 
are in reasonable agreement with experiment. Most likely the agreement with the 
p(775) could be improved upon by using a parameterization of v(z) which rises more 
rapidly to its asymptotic value at large z. Nevertheless, the purpose of this paper 
is to incorporate QCD-like chiral symmetry breaking in soft-wall AdS/QCD models, 
not just to fit data. 

3.3 Axial- vectors 

Unlike the vector field, the axial-vector couples to the scalar field VEV, producing a 
z-dependent mass term in its equation of motion. 



[19]. 



3.2 Vectors 
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Figure 5: Comparison of the predicted vector mass eigenvalues using the tanh form 
(9) ofv(z) (solid) with the QCD p mass spectrum [19]. 



The expression (16) for the axial- vector field matches that of the vector field except for 
the additional term, g 2 v 2 (z)L 2 /z 2 . Using the boundary conditions linx^^o o n (^o) — 0, 
d z a n (z — > oo) = 0, the axial-vector meson spectrum is obtained for the fixed values 
of A, m q , S, and k. 

The limiting behavior of v (z) as z — > oo leads to a constant shift between the 
vector and axial-vector spectra at high mass values. Comparing the equations of 
motion (15) and (16) for these fields one finds the asymptotic behavior 

Am 2 = {m 2 An - <)_ = g 2 ^v\z - oo) = (17) 

Together with the slope A, this determines the numerical value to be k ~ 30, although 
the best visual global fit to all the data suggests k = 15, which is the value used here. 
This is probably due to the small number of radial excitations to which we are fitting. 
The results of our analysis are plotted in Fig. 6. The ai(1260) resonance is predicted 
to within 5% and there is good agreement with the higher resonances of a±. 

Note that from (17), Am 2 > implies that k > 0, which means that the poten- 
tial in (2) is unbounded from below. To address the stability of the gravity-dilaton 
background requires a complete fluctuation analysis generalizing the work in Ref. 
[22]. Even though this is beyond the scope of the present work it does suggest that 
higher-order terms will be needed for stability. 
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Figure 6: Comparison of the numerical results for the axial-vector mass eigenvalues 
using the tanh form (9) ofv(z) (solid) with the QCD a\ mass spectrum [19]. 



4 Pion Coupling 

To confirm that our model is consistent with AdS/QCD model predictions, we calcu- 
late the pion decay constant using the formula [5] 



2 f cUo(0, 



9i ~° z 
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Here Ao(q,z) is the axial- vector bulk-to-boundary propagator with boundary condi- 
tions Aq(0, e) = 1 and d z Ao(0, z — > oo) = 0. The above formula returns the measured 
value of the pion decay constant if m q = 9.75 MeV, k = 15, and A = 0.1831 GeV 2 . 

Once we calculate the g Pn7T n, the space-like pion form factor can easily be de- 
termined from a sum over vector meson poles. However the sum converges slowly 
and numerically it is much better to use the expression in terms of the vector and 
axial-vector bulk-to-boundary propagators as in Ref. [12] 

FM 2 )= J ^e-^^^(^(^)) 2 + ^(^)-^)) 2 ), (19) 

where V (q, z) is the vector bulk-to-boundary propagator. The results of our F w (q 2 ) 
calculation are plotted in Fig. 7, and shows a slight improvement in matching the 
experimental values compared to that obtained in Ref. [12]. 



5 Conclusion 



We have shown how to incorporate chiral symmetry breaking into a soft-wall version of 
the AdS/QCD model with independent sources for spontaneous and explicit breaking. 
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Figure 7: The line shows the predicted space-like behavior of the pion form factor 
F v (q 2 ) which is compared to the experimental data obtained from [12]. The triangles 
are data from DESY, reanalyzed by [24]. The diamonds are data from Jefferson Lab 
[25]. The circles [26] as well as the star [27] are also data obtained from DESY. 



This is achieved by introducing a quartic term in the potential for the bulk scalar field 
dual to the quark bilinear operator qq. This enables us to obtain reasonable agreement 
within 10% of the QCD meson mass spectra for scalar, vector and axial- vector fields, 
although the prediction for the lowest lying p is not as good. 

Further progress must recognize the limitations of this type of phenomenological 
model. Genuine stringy behavior is most likely required to fully describe the charac- 
teristics of QCD. On the theoretical side it would be interesting to further understand 
the soft-wall model from the top-down including finding a dynamical solution of the 
features exhibited in our model along the lines considered in Ref. [15]. In addition 
the stability of the scalar potential will most likely require higher-order terms that 
can only be studied from the top-down. It is interesting that the simple 5D model 
contains QCD-like features and suggests that a further understanding of QCD can be 
obtained from the gauge/gravity correspondence. 
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